Momentum interferences of a freely expanding Bose-Einstein condensate in ID due to 

interatomic interaction change 
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A Bose-Einstein condensate may be prepared in a highly elongated harmonic trap with negligible 
interatomic interactions using a Feshbach resonance. If a strong repulsive interatomic interaction is 
switched on and the axial trap is removed to let the condensate evolve freely in the axial direction, a 
time dependent quantum interference pattern takes place in the short time (Thomas- Fermi) regime, 
in which the number of peaks of the momentum distribution increases one by one, whereas the 
spatial density barely changes. 



PACS numbers: 03.75.Kk, 03.75.Be, 39.20.+q 

The dynamics of Bose-Einstein condensates has been 
much studied, both experimentally and theoretically, for 
determining the properties of the condensate and its 
transport behavior in waveguides or free space with po- 
tential applications in nonlinear atom optics, atom chips 
and interfcrometry. Fully free (three dimensional) or di- 
mensionally constrained expansions have in particular 
been subjected to close scrutiny to obtain, from time 
series of cloud images and the appropriate theoretical 
models, information on the condensate and/or its con- 
fining potentials 1]. In reduced dimensions, the expan- 
sions have also been examined to characterize and iden- 
tify different dynamical regimes (the mean-field domi- 
nated Thomas- Fermi limit, quasi-condensates and the di- 
lute Tonks-Girardeau gas of impenetrable bosons 0,0), 
or to investigate statistical behavior in apartial release of 
the condensate into a box of finite size |j] . An effectively 
one dimensional (ID) Bose gas may be realized experi- 
mentally by a tight confinement of the atomic cloud in 
two (radial) dimensions and a weak confinement in the 
axial direction so that radial motion is constrained to 
the ground transversal state. Expansions are quite gen- 
erally described and observed in coordinate space but 
very interesting phenomena occur in momentum space 
0, Also interferences, which show the wave nature 
of the condensate and may form the basis of metro- 
logical applications, are usually apparent in coordinate 
space, between independent condensates or as a self- 
interference but, as in our present case, they may 
be genuinely momentum space effects, possibly with an 
indirect and less obvious spatial manifestation. (Other 
striking example of quantum momentum-space interfer- 
ence phenomenon for an ordinary, non-condensate, one- 
particle wavefunction colliding with a barrier has been 
discussed recently 0.) The expansions may be manipu- 
lated in different ways, e.g. by controlling the time de- 
pendence and shape of the external trapping potentials 
or by varying the interatomic interaction using a mag- 



netically tunable Feshbach resonance ^3] or an optically 
induced Feshbach resonance [Tl| . 

In this paper we shall take advantage of these con- 
trol possibilities and study a quantum interference ef- 
fect originating from a change of the interatomic interac- 
tion strength. The preparation of the condensate may be 
carried out in a ID harmonic trap with negligible inter- 
atomic interaction strength such that the ground state 
is approximately Gaussian. Removing the trap under 
these conditions would preserve the momentum distri- 
bution, but if the interatomic interaction is immediately 
increased as the trap is removed, the momentum distri- 
bution evolution changes dramatically and in a highly 
non-classical way: a time dependent interference occurs 
consisting of an orderly, one-by-one increase of the num- 
ber of peaks in the momentum distribution, see Fig. ^ 
during the early stage of the ID expansion, which is de- 
scribed accurately by the Thomas-Fermi (TF) model. In 
the same period of time of the snapshots shown, the 
spatial profile has barely evolved from the initial pro- 
file. Classical mechanics only explains the global broad- 
ening of the momentum distribution due to the release 
of (mean field) interatomic potential energy, notice the 
motion of the outer peaks in Fig. but not the oscilla- 
tory pattern, which will require a quantum interference 
analysis. A way to make the momentum interference vis- 
ible in coordinate space is to switch off the interatomic 
interaction again so that the subsequent free flight maps 
the momentum peaks into spatial ones. 

Let us now discuss the details. Assume that an effec- 
tively ID Bose-Einstein condensate is prepared in a har- 
monic trap. The condensate wave function is the ground 
state of the ID (stationary) Gross-Pitaevskii equation. 
We shall assume first that the initial interatomic inter- 
action is zero (this will be relaxed later on). Then the 
Gross-Pitaevskii equation becomes a linear Schrodinger 
equation so that the ground state condensate wave func- 
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FIG. 1: Wave function in momentum space; exact result 

i ~ 1 2 i ~ 1 2 

\ip(t, v)\ (lines), TF approximation \ipTF{t, v)\ (crosses), 
classical result P(t,v) (dashed lines); (a) t — 0, (b) t — 0.2 ms, 
(c) t = 0.4 ms, (d) t = 0.6 ms. 



tion ipo( x ) is a Gaussian, namely 



tp (x) 



muj x 



■ exp 



in 



where oj x is the axial (angular) frequency. For the rest 
of the paper we choose m = mass(Na) and oj x — 5/s. At 
time t = the confining trap is switched off and the inter- 
atomic interaction is switched on immediately. (A more 
realistic switching procedure needing a finite time will be 
discussed below.) The time evolution is now described by 
the time-dependent Gross-Pitaevskii equation, 



if i 



dtjj(t, x) 



U d 2 ib(t, x) h , , . . ,o , , , 



(2) 

where g\ is the effective ID coupling parameter. We 
are examining short times t < t c := 1/(200 0^) = 1ms 
(the factor 1/200 is arbitrary, it should be true that 
t l/w x ) for which the absolute square of the Gaus- 
sian in Eq. Q under the free evolution (i.e. Eq. (J2J 
with gi = 0) is nearly not changing in coordinate space. 
Even with g\ = 234.4 cm/s and for times t < 1ms, 
\ip{t, x)\ is nearly not changing. In contrast, the mo- 




FIG. 2: Velocities of the peaks of the wave function | ^(t, v) | : 
go = (solid lines), g — 0.002344 cm/s (dashed lines), go = 

I - I 2 

0.02344cm/s (dotted lines); TF approximation U/) rj r(t,v) 
(circles). In this and similar figures we always apply a thresh- 
old level such that the wave function is assumed to be zero if 
\tp{x)\ 2 < 0.001 * max,- IV^')! 2 - 



mentum space distribution changes substantially: more 
and more peaks are created as time increases, see Fig. ^ 
with tp(t, v) = J dx ip(t, x) exp {-i^fx) . The peak 

creation can be also seen in Fig. |3 (solid lines), where the 
velocities of the peaks versus time are plotted forming a 
characteristic structure where bifurcations and creation 
of a new central peak follow each other. 

To understand this effect, we shall approximate Eq. 
©. If we neglect the kinetic energy (TF regime) we get 

ih—ipTF{t,x) = -gi \ip T F(t,x)\ 2 ip TF (t,x). (3) 



The solution is ipTF(t,x) = V^^exp ^— \tg\ \\\)q{x)\ 
or, in momentum space, 



muj x 



dx exp 



muj x 2 



i jmuj x ( muj x 2 \ vm 
-- 9i m—ex P (-—x)- i —o 



(4) 



Note that in the TF regime the spatial density remains 
unchanged. The results for Eq. |0J are also plotted in 
Figs, ^and |21 There is a good agreement between the 
exact result and TF, and both show the same interference 
behavior. It is clear that the nonlincarity is playing a key 
role in the effect. 

We may compare the quantum dynamics with a similar 
classical one. Let us assume an ensemble of classical par- 
ticles where the probability density of initial positions 

and velocities is given by po(x,v) = |i/>o(:r)| 2 ^o(^) 

with ipo{v) given by Eq. Q). The probability den- 
sity is then evolved with the classical Liouville equation 
with the Hamiltonian H — ^gi J dv' pt(x,v'), without 
kinetic term in analogy to the quantum Hamiltonian of 
the TF regime, see Eq. ©. The solution is given by 
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FIG. 3: Wave function for v = 0: exact result \^f(t, 0)| 

I " 1 2 

(lines), TF approximation \&tf(£, 0) (crosses), stationary- 

l " l 2 

phase approximation w s (i) 0) (dashed line); the filled circle 
marks to. 



FIG. 4: Velocities of the peaks of the exact result $ 

I " 1 2 

(lines), the stationary- phase approximation \ips\ (triangles); 
the outer thick-dotted lines mark the critical values ±v c (t). 



pt(x,v) = po(x,v + at) with a = ffj^ Jdvp {x,v). 
P(t,v) :— f dxp t (x,v) is also plotted in Fig. ^ The 
classical picture provides, approximately, the outer peaks 
moving outwards with time because of the release of po- 
tential energy, but not the central ones, which are there- 
fore associated with quantum interference. Let us exam- 
ine this interference. 

Defining dimcnsionless quantities ( = \J mLu x /h x, t = 
\J muj x /h g\t, and k — yfm/%LO x v, we write Eq. @ as 



^tf{t, k) = 



d£ exp 



1 



2 



muj x 



V20F 



< 2 + -c 



(5) 



An important characteristic time is the value To when the 

2 



second maximum of 



in the TF approx- 



^tf(to, k = 0) 
imation appears (the first maximum is at r = 0), see Fig. 
13 This sets the scale of the oscillations and the value is 
found numerically, 



t fa 27.703 



to 



27.703 

91 



This time t Q should be much smaller than the critical time 
t c when the kinetic energy starts to influence. Therefore 
we get a necessary condition for g%, 



tc 



1 



200 uj x 



> in 



27.703 

TlLUr, 



mu) x 



gx > 27.703 x 200W — - pa 65.12 cm/s 



m 



We will simplify further Eq. (J5J by means of the station- 
ary phase method. If k/t is constant and t — * oo then 
the main contribution to the integral @ comes from the 
values ( which fulfill 



d_ 



2^ 



-c 

r 



= o^-c- 



1 = o. 

T 



If < \h/t\ < exp (— 1/2) /v2tt there are two solu- 
tions Co, Ci witn < | Co | < 1/V2 < |Ci|- Physically 
these are two positions in which the force exerted on the 
atom is equal because the Gaussian profile of the po- 
tential changes from concave-down around the center, to 
concave- up in the tails. Thus these two positions con- 
tribute to the same momentum and the corresponding 
amplitudes will interfere quantum mechanically. The sta- 
tionary phase approximation of Eq. (J2J) is 



rriuj-. 



1 



(6) 



exp 



exp 



resulting from 



Fig. 0] shows also the peaks of 

using Eq. © in the range < \k/t\ < exp (—1/2) /V^tt. 
Except for the outer peaks and the k = line, Eq. © 
gives the correct peak behavior. The opening of the cone 
of the effect can be approximated by the definition of a 
critical k c (t) or v c (t) to make the interference possible, 



Kc _ CXp (-1/2) 



exp (-1/2) 

K c {t) = — T, 



2tt 



or v c (t) = exp (— l/2)u) x git/y/2TT. 

The interference pattern in momentum space can be 

seen in coordinate space if the interaction is switched off 

2 

at a given time i //- Then ?p(v,t) is not changing and 
therefore the position and number of peaks is not chang- 
ing. On the other hand, the different peaks separate in 
coordinate space. Therefore for different times t f / a dif- 
ferent peak pattern will appear in coordinate space (see 
Fig. 03). Only the central peak cannot be seen. 

Up to now we have considered an abrupt switching of 
the interaction, but the effect remains for switching times 
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FIG. 5: Wave function in coordinate space; t = 40 ms; At = 
(solid lines), At = 0.1ms (dashed lines); the positions of the 
maxima are marked by crosses for At — 0; (a) t a ff = 0.2 ms, 
(b) t ff = 0.4 ms, (c) t ff = 0.6 ms. 



At even of the order of to . For a time dependent coupling 
parameter given by 



g(t)=gix 



/(*) 

f(toff - t) 



t<0oit> t of f 

At <t< t off - At 

< t < At 

toff — At < t < t ff 



(7) 



with /(<) = {t/At) 2 (3-2t/At), the result for At = 0.1 ms 
is shown in Fig. [SJ the effect is not changing qualitatively, 
only the positions of the peaks are squeezed. 

We shall also examine the stability with respect to a 
non-zero value of the coupling constant used to prepare 
the initial state, go- The initial ground state is then no 
longer a Gaussian and can only be calculated numerically. 
The effect survives as long as 50 "C g\ but the interference 
pattern is again squeezed, see Fig. |21 

Summarizing we have examined the short-time behav- 
ior of the evolution of a Bose-Einstein condensate in 
ID when the interatomic interaction is negligible for the 
preparation in the harmonic trap, and strongly increased 
when the potential trapping is removed. We have found 
a quantum interference effect in momentum space origi- 



nated from the interatomic interaction change. The mo- 
mentum distribution expands due to the release of mean 
field energy and the number of peaks increases with time 
because of the interference of two positions in coordinate 
space contributing to the same momentum or velocity. 
The effect is stable in a parameter range and and could 
be observed with current technology. 
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